Classical heisenberg antiferromagnet away from the pyrochlore lattice limit: 

entropic versus energetic selection 
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The stability of the disordered ground state of the classical Heisenberg pyrochlore antiferromagnet 
is studied within extensive Monte Carlo simulations by introducing an additional exchange interac- 
tion J' that interpolates between the pyrochlore lattice (J' = 0) and the face-centered cubic lattice 
(J' = J). It is found that for J'/J as low as J' /J > 0.01, the system is long range ordered : the 
disordered ground state of the pyrochlore antiferromagnet is unstable when introducing very small 
deviations from the pure J' — limit. Furthermore, it is found that the selected phase is a collinear 
state energetically greater than the incommensurate phase suggested by a mean field analysis. To 
our knowledge this is the first example where entropic selection prevails over the energetic one. 

PACS numbers: 75.10Hk, 75.40Cx, 75.40 Mg 



I. INTRODUCTION 

In recent years, there has been a great interest in 
the study of fully frustrated Heisenberg antiferromag- 
nets (HAF), both from the experimental and theoretical 
point of vievcrEI. These systems often show unconven- 
tional groiind-states, noncollinear, incommensurate^ or 
disordered!! In these systems, the classical Mean Field 
(MF) description indicates a macroscopic ground-state 
degeneracy and no long range order is predicted. This 
property is directly related to the extensive number of de- 
grees of freedom that still fluctuate at low temperatureB. 
Nevertheless, in many cases, thermal or quantum fluctu- 
ations give rise to an entropic selection within-the ground 
state manifold, known as "order by disorder"Q. This phe- 
nomenon may induce a selection of a subset of the contin- 
uous degenerate-ground states as it does for the kagome 
antiferromagnetD, and sometimes can even select a par- 
ticular ground state as it iSpthe case for the face centered 
cubic (fee) antiferromagnettl 

Among these systems, the 3D pyrochlore lattice is of 
particular interest since it does not display such "or- 
der by disorjjLe£' and remains disordered down to zero 
temperatureU'tm. It consists in a three dimensional ar- 
rangement of corner-sharing tetrahedra (see Fig. Q). All 
magnetic compounds which cristallizc in the pyrochlore 
structure exhibit unusual magnetic properties. A few 
order at well defined Neel temperaturea while many of 
them behave as. spin glasses despite the high degree 
of stoichiometryu. In the Heisenberg antiferromagnetic 
5 = 1/2 quantum case, the pyrochlore lattice exhibits 
no magnetic transition even at T=0K, and behaves as 
a spin-liquid. The spin-spin correlation length never ex- 
ceeds the inter atomic distance and it is expected that 
there is a cantinuum of singlet states within the singlet- 
triplet gapO. For classical Heisenberg spins, Monte Carlo 




FIG. 1: The pyrochlore lattice. 



simulations have shown that the pyrochlore HAF also be- 
haves as a spin-liquid, with short ranged spin correlations 
both in space and timeOiJ, which allows for its description 
to take into account only the tetrahedral unit cellid. The 
exactly solvable infinite-component spin vector model on 
the pyrochlcpe lattice also does not order down to zero 
temperatureta. Classical models providing much infor- 
mations close to the characteristics obtained in the quan- 
tum case, it indicates that the low-energy physics in the 
pyrochlore structure is mainly related to the geometry of 
the lattice. 

An interesting aspect of the problem is to study 
how this disordered classical ground-state resists to per- 
turbations as it could enlarge its observation in real 
compounds. It has been shown that ordering is very 
easily induced, either by including furttuiE- neighbour 
interactionliiliJ or easy plane anisotropyElE£l or lattice 
deformationala. The purpose of this paper is to explore 
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the effect of coupling the spins of the pyrochlore net to 
the remaining nearest-neighbor (n.n) spins in the fee net 
via standard classical Monte Carlo (MC) calculations. 
The question is wether the system will select an ordered 
state when switching on this coupling or maintain disor- 
der as in the pyrochlore HAF. As obtained in previous 
works, an energetic comparison of all magnetic phases at 
the MF level already lifts the degeneracy and strongly 
supports that any deviations from the pyrochlore limit 
will drive the system to an ordered state (sec. |l]). But 
strikingly, it is shown that thermal fluctuations do not 
select the expected mean field ground state but prefer a 
collinear phase which is more stable from the entropic 
point of view (sec. III). 



II. FROM THE PYROCHLORE LATTICE TO 
THE FCC LATTICE : THE J - J' MODEL 

In both fee and pyrochlore lattices, the classical en- 
ergy is minimized by any configuration for which the to- 
tal spin of each elementary tetrahedron in the lattice is 
zero. The main difference in the ground-states of these 
two lattices comes from the fact that in the fee net the 
ordering at one tetrahedron will select a unique global 
ordering, while in the pyrochlore net, since the tetrahe- 
dra are more sparsely connected, there remains an exten- 
sive number of degrees of freedomcl. In the fee HAF, the 
MF ground-states are degenerate along the ir/a(l,q, 0) 
direction in reciprocal space. Thermal fluctuations break 
this degeneracyp.and select the collinear ordered phase 
qo = (ir/a, 0, 0)H. These results have been confirmed by 
standard Monte Carlo simulationalU which show evidence 
for a first order phase transition. In the pure pyrochlore 
HAF, the MF ground-states are degenerate throughout 
the entire Brillouin zone. Numerical works have shown 
that no long-range order occurs at nonzero termperature 
in that case, even in a regime where the .-temperature is 
below the energy scale set by interactionBEI. In between 
these two limits, we have considered the J — J' fee HAF, 
where J couples the n.n. spins on a pyrochlore net and 
J' the remainder n.n. spins on the underlying fee net (see 
Fig. ||). This leads to the following classical Heisenberg 
hamiltonian : 
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where Si is a vector spin of unit length occupying the 
i th fee lattice site, i and i' design sites on two pyrochlore 
nets shifted by (2a, 0,0) respectively and the sums <> 
run over n.n. pairs. We have limited our study to the 
antiferromagnetic case, i.e. all negative interactions. J — 
J' < corresponds to the fee HAF and J' = 0, J < 
leads to the pyrochlore HAF. The cubic lattice constant 
is 2a. 

We have first studied the MF ground-states degeneracy 
for small a = J' / J. One way to do so is to compute 
the Fourier transform of the interactions on the lattice, 




FIG. 2: The J — J' n.n. fee lattice. The fee lattice is decom- 
posed into two pyrochlore nets : A sites (plain lines) and B 
sites (dashed lines). J couples A — A n.n. sites (plain lines) 
while J' couples B — B n.n. sites (dashed lines) and A — B 
n.n. sites (not drawn for clarity). This model contains 8 sites 
per unit cell. 
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FIG. 3: Largest eigenvalue A maa; (q)/|J| along the high sym- 
metry axis of the fee Brillouin zone for a — 0.01 < a c ~ 0.21. 



,/ q . The ground-state (if unique) then corresponds to the 
highest eigenvalue A mQ piqi ) of J q and is described by the 
propagation vector qiE3. For the J — J' model, J q is a 
8x8 matrix that has already been calculalted by Reimers 
et alH. 

For < a < a c ~ 0.21, we have reported the largest 
eigenvalue A maa; (q)/|J| along the high symmetry axis of 
the fee Brillouin zone on figure ^. The degeneracy is com- 
pletely lifted and A maa; (q) has an absolute maximum qi 
along the (0, q, q) direction, selecting an incommensurate 
phase of wave vector qi . Within the MF approximation 
the J— J' model selects an ordered phase for < a < a c . 

For a > a c , the situation is very similar to the fee 
case. The MF degeneracy is lifted except on lines along 
7r/a(I,g,0) direction : the system remains disordered in 
that case. 

We have reported in figure | the MF energy difference 
AS/) J| between the two ground states: the qi incom- 
mensurate state and the 7r/a(l, q, 0) states (including the 
qo collinear state) as a funtion of a. We see that AE/\ J 
is much lower than the energy scale set by the interac- 
tion (AE/\J\ ~ a/100). Thus the degeneracy is very 
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TABLE I: Extrapolated zero temperature specific heat for 
various values of a = J' /J from a = (pyrochlore HAF) 
to a = 1 (fee HAF). This system displays long range 
order(C = 1) for all a > 0. 



0.4 0.6 . 



a = J' /J 


C in kB units 


(pyrochlore HAF) 


0.73±0.03 


0.01 


1.01±0.02 


0.1 


1.00±0.02 


1 (fee HAF) 


0.998±0.005 



FIG. 4: MF energy difference 10 3 x AE/\J\ between the two 
MF ground-states as a funtion of a (AE — E(qi) — E(qo), 
see text). 



slighty lifted for a < a c : the highest branch is almost 
flat over the whole Brillouin zone and very similar to the 
pyrochlore one. If we add thermal fluctuations on these 
MF ground-states, it is not clear what will happen. They 
can either select a particular state, the qi noncollincar 
state (energetic selection) or a collinear state (entropic 
selection) or induce disorder. 



III. MONTE CARLO RESULTS 

We have performed standard MC simulations using the 
sample sizes of N = 864—6912 spins with periodic bound- 
ary conditions, over a temperature range extending down 
to T/\J\ — 10~ 4 . We started from both random and 
ordered initial configurations. In each run, up to 10 7 
MC steps per spin have been discarded for equilibrating 
and up to 2 x 10 7 MC steps per spins have been used 
for averaging. We have checked that our simulations n 
produce the Monte-Carfor*esults for the pyrochlore netl 
and those for the foe netO. In order to determine un- 
ambigously whether the J — J' foe HAF order or not, we 
have measured more specifically the specific heat per spin 
and a correlation function which quantify the collinearity 
of spins. 

An unbiaised way to look for fluctuations induced or- 
der is to measure heat capacity, since it is very sensi- 
tive to the presence of zero modes (modes whose energy 
is independent-jOf displacement to second order at all 
wave vectors)l§El. From the classical equipartition theory 
quadratic modes contribute a factor T 1 / 2 to the parti- 
tion function Z, while quartic modes contibute a factor 
T 1 / 4 : at low temperature, the behavior is thus domi- 
nated by fluctuations around the ground-states with the 
largest number of zero modes. These zero modes leave a 
signature in the heat capacity (since they contribute to 
&b/4 instead of fce/2 for the quadratic modes) that can 
be checked with Monte Carlo simulations. 

In the n.n. pyrochlore HAF, fluctuations around the 
collinear states lead to one zero mode per tetrahedron re- 
sulting in a heat capacity per spin C = 3fce/4 in absence 



of order (instead of C = 5fcs/8 if a collinear state is se- 
lected). This value-has already been checked by Monte 
Carlo calculationsclu. In the fee n.n. HAF, fluctuations 
select the collinear state and since it contains only degen- 
eration lines, the heat capacity per spin is simply equal 
to C = in that case. 

As the heat capacity derived from energy fluctuations 
by MC simulations is rather noisy at low temperatures, 
we have estimated the heat capacity directly from the 
slope of the asymptotic internal energy at low tempera- 
tures. The obtained results are reported in Table | and 
clearly show that for a > 0.01, thermal fluctuations se- 
lect an ordered state at zero temperature. To show more 
directly whether the system selects the collinear ordering 
or not, we have calculated a parameter that measures the 
collinearitytj : 



P(n,j) = 7; ( < (s,.s,) 2 > 



(2) 



where r^j is the distance between i and j sites. So P=l 
for collinear spins and P=0 in the high temperature limit. 
This parameter is represented in figure ^| as a function of 
the separation r (in units of n.n. distances) at T/\J\ = 
5 x 10~ 4 for different values of a. We find that as soon 
as a > 0.01, the collinearity parameter is of long range, 
and close to 1, contrary to the pyrochlore limit which 
shows a correlation for r < 2 only. We can note that 
the order parameter P(r — > oo) is rather lower than 1 
despite the low-temperature takers as for the XY model 
on the pyrochlore antiferromagneifl. Figure ^ shows the 
temperature dependence of the collinearity parameter for 
neighbouring spins r = 1. As soon as a > 0.01, the 
spins become collinear in the low temperature limit, by 
contrast with the pyrochlore case which has a small but 
finite correlation length. We have reported in figure ^ the 
critical temperature as a function of | J'| for J' ^ 0. It is 
of the order of | J'|/2. We have checked that these results 
do not depend on the initial configuration choosen, by 
starting from collinear and random states. And the sizes 
considered are large enough since no marked finite-size 
effects arc observed. 

We have checked the low temperature phase starting 
from collinear and random states. We found the qo 
collinear state for a > 0.01. Let's notice that starting 
from a random configuration, stacking faults set in at 
low temperatures in MC calculations, as it is usual in 
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FIG. 5: Dependence of the collinearity parameter P(r) at 
T/\J\ = 5 x 1(T 4 for different values of a = J'/J, 4 x L 3 , 
L — 12 Heisenberg spins and starting from random spin con- 
figurations, r is given in units of n.n. distances. 

fully frustrated systems^. At T/\J\ = 5 x 1CT 4 and for 
a = 0.01 and N = 6912, we found two collinear blocks 
(of qo type) of around N/2 spins and two non collinear 
walls of few hundred spins. For a > a c , we observe an 
order by disorder effect very similar to the fee HAF : 
thermal fluctuations induces the collinear magnetic or- 
der qo via an entropic selection. For < a < a c , the 
low temperature state selected by thermal fluctuations 
is not the qi incommensurate state selected at the MF 
level, neither a state with a MF energy close to the qi 
state and commensurate with the periodic boudary con- 
ditions we have taken, but the collinear state qo whose 
wave vector is far from qi in the Brillouin zone. Even for 
a « 0.1, where the MF energy difference between these 
two ground states becomes greater (AE/\J\ ~ 10~ 3 , see 
Fig. 0) than the lowest temperature we have considered 
T/\J\ = 5 x 10~ 4 , the system still selects the collinear 
state. The entropic selection always overcomes the ener- 
getic one in that case. 



one. This is the first case, to our knowledge, where such 
a phenomenon is pointed out. The pyrochlore lattice 
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FIG. 6: Temperature dependence of the collinearity parame- 
ter P(r = 1) for different values of a = J'/J, different sizes 
4 x L 3 and starting from random spin configurations. 
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FIG. 7: Phase diagram of the J — J' model versus a — J'/J 
(qo = (7r/a,0,0)). 



IV. CONCLUSION 

We have studied in this paper the effect of small de- 
viations from the pyrochlore n.n. HAF to the fee one 
via standard MC simulations. As soon as we switch on 
J' (a = J'/J > 0.01), correlation functions are long- 
ranged in space and the specific heat at zero tempera- 
ture does not indicate the presence of zero modes. The 
J — J' HAF selects the qo = (tv/ci, 0, 0) collinear state at 
low-temperature and ordering occurs at finite tempera- 
ture. For small a, thermal fluctuations do not select the 
MF ground-states at low temperatures but the collinear 
state. The entropic selection prevails over the energetic 



thus remains the only 3D example known with classical 
Heisenberg spins remaining disordered at all tempera- 
tures. 
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